A new universality class describing the insulating regime of disordered wires with 

strong spin-orbit scattering 
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We discuss the distribution of transmission eigenvalues in the strongly localized regime in the pres- 
ence of both a magnetic field and spin-orbit scattering. We show that, under suitable conditions, 
this distribution can be described by a new universality class labelled not only by the index but 
also by a new index r\. This result is obtained by mapping the problem into that of a suitable 
Calogero-Sutherland model. 
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Magnetoconductance (MC) in the strongly localized 
regime of mesoscopic systems has been the subject of an 
increasing interest during last years jjj - ||. While the 
same phenomenon in the weakly localized regime is by 
now rather well understood, and a good agreement has 
been reached among different theoretical descriptions and 
between their predictions and the experimental results, 
the situation is less clear in the strongly localized regime. 
In particular it is not clear to which extent the Random 
Matrix Theory (RMT) of quantum transport and 
the corresponding classification in terms of universality 
classes, originally developed to describe the weakly lo- 
calized regime, can also be extended to the insulating- 
regime. Indeed, it seems that two different insulating 
regimes exist. The first one (strongly insulating regime) 
in which the localization length £ (to be defined below) 
is much smaller than the hopping length th , seems to 
be well described by a directed-path approach (see ref. 

and references therein). The second one (weakly 
insulating regime), close to the Anderson transition, in 
which th/£ ~ 1 and backscattering paths are important, 
seems to be well described within a RMT framework. 
This would be an important result due to the simplicity 
and generality of the RMT approach: the input required 
is only the symmetry of the Hamiltonian and the result- 
ing classification is completely encoded by a single index 
0. The above scenario is however complicated by the 
presence of some conflicting theoretical analyses. The 
disagreement between the various approaches becomes 
evident when looking at the MC behaviour in presence 
of strong spin orbit scattering (SOS). In fact in the in- 
sulating regime the localization length is related to the 
conductance by: g = exp(— 2L/£), (L=length of the sam- 
ple). In the standard RMT approach £ is proportional 
to and, as a consequence, the magnetoconductance is 
directly related to the the change of as a function of the 
magnetic field. The RMT approach based on the maxi- 
mal entropy ansatz predicts that the effect of a magnetic 
field on a sample with SOS is to change from = 4 
( "symplectic ensamble" : spin-rotation symmetry broken 
but time reversal symmetry conserved) to = 2 ("uni- 



tary ensamble" : time reversal symmetry broken, status of 
the spin-rotation symmetry irrelevant). Defining as £{H) 
the localization length in presence of a magnetic field, 
from the above analysis it follows that£(ii")/£ = 1/2, and 
we expect a strong negative MC j|] . On the contrary an 
exact calculation within the non-linear a model [Q shows 
that if the spin-rotation symmetry is broken, remains 
equal to 4 even in presence of a breaking of the time re- 
versal symmetry (thus implying a zero MC) and that in 
this case the ensamble behaves as unitary (0 = 2) only 
with respect to the level statistics. In this letter we shall 
try to clarify this issue. Let us state our main results: 

a] The main point we want to stress is that, while 
in the weakly localized regime all the relevant physical 
properties are completely determined by the level statis- 
tics (namely the value of 0), when extending the RMT 
approach to the insulating regime, due to the fact that 
in this case the conductance is dominated by the lowest 
eigenvalue all the details of the chosen RMT model (num- 
ber of degrees of freedom, possibly the presence of new 
critical indices) become important. 

b] In particular this gives a way to reconciliate the 
above described controversy. We shall show below that 
also within the maximal entropy approach a zero MC 
should be expected in the insulating regime. 

c] It is however very interesting to notice that in 
some recent experiments [Q on samples characterized by 
strong SOS a clear negative MC has been observed in the 
weakly insulating regime, in contradiction with the above 
results. We shall show below that under suitable hypoth- 
esis, it does exist a new universality class in the RMT 
description of quantum transport which indeed admits 
negative MC in the insulating regime. This new class 
has = 4, as in S, but a new index rj exists that allows 
to distinguish this case from the ordinary symplectic en- 
samble. As a consequence of this new 77 index a negative 
MC is expected even if = 4. 

To show these results let us first briefly remind few 
known features of the RMT description of quantum 
transport. The central object in this approach is the 
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2N x 2N Transfer Matrix M (where N is the number 
of scattering channels at the Fermi level) which relates 
the incoming and outgoing fluxes in a two-probe geom- 
etry (namely a finite disordered section of length L and 
transverse width L t , to which current is supplied by two 
semi- infinite ordered leads): 

"(£)-(?) • M 

where the 1,0,1', O' vectors describe the incoming and 
outgoing wave amplitudes on the left and right respec- 
tively. To construct the RMT one must first impose the 
physical requirements which define the ensamble of the 
matrices: for instance flux conservation, or time conjuga- 
tion, if time reversal symmetry is conserved. This defines 
the group G to which the random matrix belongs. Then 
one assumes that the physics of the problem only de- 
pends on a subset of the degrees of freedom of the matrix. 
It turns out that a good parametrization in our case is 
given by the N real, positive eigenvalues A^ of the matrix 
Q = \ [M f M + (M^M)' 1 - 2] . These parameters have an 
important physical interpretation since they are related 
to the transmission eigenvalues Tj by A; = (1 — T,-)/Tj in 
terms of which all the interesting physical properties can 
be evaluated. For instance the conductance g is given 
by 9 = 9o^2 n T n (with go — 2e 2 /h). This parametriza- 
tion allows to identify a subgroup H of invariance in our 
ensamble. The matrix ensamble is thus related to a suit- 
able G/H coset. It turns out that in all the relevant 
cases G/H is actually a symmetric space. The last step 
is the identification of a dynamical principle to describe 
the "time" evolution of the eigenvalues. This will allow to 
construct explicitly the probability distribution P({Ai}) 
for the eigenvalues A^. In the case of RMT for quantum 
transport the "time" s = L/l is actually the length of the 
conductor L measured in unities of the mean free path I; 
the dynamical principle is the so called "local maximum 
entropy principle" which simply implies the free motion 
of the eigenvalues in the symmetric space G/H and, as 
a consequence, the evolution operator is the radial part 
of the Laplace-Beltrami operator on G/H. The corre- 
sponding (coupled) differential equation for the probabil- 
ity distribution P(X) of the eigenvalues is known as the 
Dorokhov-Mello-Pereyra-Kumar (DMPK) equation || 

¥■ = - £ T^ 1 + A ») J Wx- P = D P > ^ 
os 7 f—' d\i dXi 

where 7 = 0N + 2-/3 and j3 € {1, 2, 4} is the symmetry 
index discussed above. </({A„}) denotes the Jacobian 
from the matrix to the eigenvalue space: 

j({A4)=ni^-^i /3 - ( 3 ) 

i<j 

We have seen that the RMT is completely identified by 
its G/H symmetric space. All the irreducible symmetric 



spaces of classical type can be classified with essentially 
the same techniques used for the Lie algebras. They fall 
into 11 classes labelled by the type of root system and by 
the multiplicities of the various roots |fTo|| . The DMPK 
operator D is completely fixed once the information on 
the type of root system and the root multiplicities are 
given. The Dyson ensambles, which first appeared in the 
RMT of nuclear energy levels, (and are also relevant in 
the RMT of quantum transport based on the scattering 
matrix [jll]]) are related to the symmetric spaces labelled 
by the Dynkin diagrams of type Am- There are exactly 
three classes of A m spaces, which are distinguished by 
the index (3 which labels the multiplicity of the roots. 
can take the values 1,2 or 4 and exactly coincides with 
the symmetry index of the corresponding Random Ma- 
trix ensamble. Instead, the ensambles or random trans- 
fer matrices in which we are interested coincide with the 
symmetric spaces labelled by Dynkin diagrams of type 
Cjv (see Tab.l). It is exactly the difference between Am 
and Cm type of symmetric spaces the origin of the differ- 
ences between RMT based on transfer matrices and the 
Wigner-Dyson ones. 



G 


H 


d 





V 


Sp(2N, R) 
SU(N,N) 
SO*{4N) 
USp(2N, 2N) 
S P (2N,C) 


U(N) 

SU(N)®SU(N)®U(1) 
U(2N) 
USp(2N) ® USp(2N) 
USp(2N) 


N(N + 1) 

2N 2 
N(N - 1) 

4N 2 
N(2N + 1) 


1 

2 
4 
4 
2 


1 
1 
1 

3 
2 



TABLE I. Irreducible symmetric spaces of type Cjv. In 
the first two columns the group G and the maximal subgroup 
H which define the symmetric space. In the third column 
the dimensions d of the space. All these spaces are labelled 
by the Dynkin diagram Cjv. In the fourth column the multi- 
plicity f3 of the ordinary roots of Cn- In the last column the 
multiplicity r\ of the long root. 

The most relevant feature of the Cm symmetric spaces 
with respect of the Am ones is that they are described by 
two critical indices instead of one. This is due to the fact 
that the Cm Dynkin diagrams have two types of roots, 
each with its particular multiplicity. We shall call this 
second index 77 in the following. This peculiar feature of 
these ensambles was up to now ignored because in the 
weakly localized regime it is only the value of f3 which 
matters, and also because the three cases which have 
been studied up to now (the first three lines of tab.l) 
have the same value of the second index 77 = 1, thus 
leading to the same DMPK equation (with no explicit 
1] dependence). However we shall show now that in the 
insulating regime, if both the spin-rotation and the time 
reversal symmetry are broken, also the fourth ensamble 
of those listed in tab.l could become relevant, and in this 
case the value of 77 changes. 
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As a preliminary step to prove this assertion, let us 
remind the group properties of the symplectic ensamble 
of transfer matrices (spin-rotation broken, time reversal 
conserved). In this case each one of the /, I', O, O' vectors 
is a collection of iV spinors (one for each channel); each 
spinor contains the two spin degrees of freedom. Hence 
M is in this case a 4N x AN complex matrix. 

Flux conservation implies 
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which can be written as: 



M f E,M = E. 



where 



E z = 



1 
1 



(5) 



(6) 



(7) 



where 1 and denote the 2N x 2iV unit and zero ma- 
trix. Eq.(||) implies that the matrix M belongs to the 
unitary group U (2N, 2N) group (since it preserves the 
metric E z ). Time reversal must then be imposed in two 
steps. At the level of the single spinor we must require 
invariance with respect to the time reversal operator: 
9 = ia y C where C is the complex conjugation opera- 
tor and a y is the Pauli matrix: 



1 

-1 



C 



(8) 



In addition, at the level of the whole matrix M one must 
impose the following constraint 



with 



M* 



Ex — 



E T ME, 



1 

1 



(9) 



(10) 



which, together with cq.lQ) implies that M € 
SO* (4 N). Let us now switch on the magnetic field. We 
have already seen that if both the constraints (|J) and (||) 
are eliminated we have M € U(2N, 2N), namely that M 
belongs in this case to the unitary ensamble (3 = 2 listed 
in tab.l, but with N -> 27V. This gives the standard RMT 
prediction: (3 changes from 4 to 2 in relation to the level 
statistics, as a consequence of the application of a mag- 
netic field. However for all those observables which are 
not only related to the level statistics, but depend more 
generally on 7, the two changes (3 : 4 — * 2, N — > 2N com- 
pensate and the model effectively behaves as if (3 = 4 as 
in the a model result. This happens in particular for the 
localization length £ in the insulating regime in which 
we are interested. In fact we shall show below that £ 
is proportional to 7. As mentioned above this seems to 
exclude a negative MC in the RMT approach, in clear 



contradiction to the experimental results. To open the 
possibility to find a negative MC with RMT, we must 
make a further hypothesis. Since the two constraint (||) 
and (||) are not on the same footing, the first being a local 
and the second a global property of the electrons, they 
need not to be simultaneously fulfilled. In particular it 
could exist a range of parameters (say, for low magnetic 
fields) in which only the second one is eliminated but the 
first one survives. Physically this would mean that lo- 
cally the Kramers degeneracy is conserved, but globally 
time- reversal is broken, with all the known consequences 
on the back-scattering paths. Let us see which are the 
consequences of this hypothesis. It is easy to see that 
eq. alone implies the conservation of the antisymmet- 
ric metric and that in this case the joint application of (|^) 
and (|) leads to M <E USp(2N, 2N). This choice of con- 
straints defines the new universality class anticipated in 
the title. Since USp{2N,2N) = U{N,N,Q) we see that 
in this case M has a very natural description in terms of 
quaternions, being the most general unitary quaternion 
matrix which preserves the metric E z . Any matrix M of 
this type can be parametrized as: 



M 



■i« 




VTTa 

'a VTTa 



u (3) 









(12) 



where the u">, (i = 1, 2, 3, 4) are N x N quaternion uni- 
tary matrices and A is a N x N real, diagonal, quater- 
nion matrix with non negative elements Ai, A2, • • • Ajv. 
This parametrization is nothing else than the explicit 

construction of the symmetric space u^Q^^l Q } > 
which (through the identification of the quaternion uni- 
tary groups as unitary symplectic groups on the complex) 
is exactly equivalent to the one listed in the fourth line 
of tab.l . Let us now generalize the DMPK equation 
to take into account this new critical index. The (gen- 
eralized) DMPK equation can be obtained as a particu- 
lar case of the general Calogero-Sutherland construction 

(see Hll). The result is: 



dP 

ds 
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where 



j({K})=n\x j -\ i fi[[\ i (i+\ i )} : 



(14) 



(15) 



i<j 



7 = -^i((3(N — 1) + 1 + 77) and the indices f3 and 77 are 
listed in tab.l. 

The DMPK equation greatly simplifies if one chooses 
a new set of variables {x n } related to the eigenvalues 
{Ai} by A n = sinh 2 x n . Then the solution can be written 
in terms of so called "zonal spherical functions" $fc(#), 
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x = {xi, • • • , xn}, k = {fci, • • • , /cjv} (which depend on 
the chosen symmetric space). 



P({x n },s) = [t(x)]* I dke~^ s MX) 



|A(fc)| 2 



F(k) (16) 



where 



£(x) = Y\ I sinh 2 xj - sinh 2 x^ JJ | sinh 2x l \% ; (17) 



irk,; 



F(fc)=n% tanh (^)(fci+i) 



\Mk)\ 2 = n 



rn<j 



r(f + ^)r(f 



(18) 



(19) 



and r denotes the Euler gamma function. 

In the insulating regime {k <C 1) we can use the asymp- 
totic expansion described in Jljj for the zonal spherical 
function. In this way the integration over k can be done 
explicitly and gives Ea: 



SO* (4AQ 
U{2N) 



USp(2N,2N) 



USp(2N)®US P (2N) and £(#)/£ = 1/6 

implying a negative MC. 

Needless to say that the negative MC measured in the 
experiments could also simply indicate a failure of the 
RMT approach outside the metallic phase. Moreover 
one must be aware of the fact that all the above results 
are constrained to the quasi one-dimensional limit which 
underlies the DMPK equation. Nevertheless, it would 
be very interesting to test the possible existence of this 
new universality class. In this respect a stringent and 
simple test is given by the ratio between the mean and 
the variance of logg. This ratio is easily extracted from 
eq.s(pUE2) and turns out to be 



var (logy) 2 
<logff> v' 



(24) 



This observable gives a direct measure of r\ and could 
play in the insulating regime the role which is played in 
the metallic one by the universal conductance fluctua- 
tions for p. We hope that our results will stimulate new 
experiments in order clarify this issue. 



2 2\ 

Xj x i ) 



P({i„}, s) — Y[ | sinh 2 Xj — sinh 2 x, | 2 [( 

x |^exp(-a; 2 7/(2s))a; ^ (sinh2x ^ ) ^, 



/2 



(20) 



Ordering the Xn's from small to large and using the 
fact that in this regime 1 <C xi < 12 < ■ ■ ■ C sjv we 
can approximate the eigenvalue distribution as follows: 



N 



(21) 



where x n = 4(r? + f3(n — 1)). In the strongly localized 
limit the conductance is dominated by the first eigenvalue 
of the transmission matrix: 

<log 5 >= -H- -logT! ~ -2xx - (22) 
£ 7 



and hence: 
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(23) 



Setting now the proper values of (3, N and r\ for 
the various regimes, and taking the large N limit, 
we see, as anticipated, that if both time reversal and 
Kramers degeneracy are broken, namely if we move from 

U(2N) ) ~^ SU(2N)®SU(2N)®U(1) tnen ^ ^ 0eS n0t cnan g e i 

hence £(H)/£ = 1 and the MC is 0. If we assume that 
Kramers degeneracy is conserved then we move from 
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